Noise induced Hopf bifurcation 
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We consider eflect of stochastic sources upon self-organization process being initiated with creation 
of the limit cycle induced by the Hopf bifurcation. General relations obtained are applied to the 
stochastic Lorenz system to show that departure from equilibrium steady state can destroy the 
limit cycle in dependence of relation between characteristic scales of temporal variation of principle 
variables. Noise induced resonance related to the limit cycle is found to appear if the fastest 
variations displays a principle variable, which is coupled with two different degrees of freedom or 
more. 
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I. INTRODUCTION 

Interplay between noise and non-linearity of dynami- 
cal systems [H is known to arrive at crucial changing in 
behavior of systems displaying noise-induced [1, 3| and 
recurrence [1, [H] phase transitions, stochastic resonance 
B 01 , noise induced pattern formation Q , and noise 
induced transport (sl. Il0l|. The constructive role of noise 
on dynamical systems includes hopping between multi- 
ple stable attractors [ll|, [13] , stabilization of the Lorenz 



attractor near the threshold of its formation |13l . Il4| and 
stabilization of resonance related to the limit cycle near 
the Hopf bifurcation 13| . Such type behavior is inherent 
in systems which involve discrete entities (for instance, in 
ecological systems individuals form population stochasti- 
cally in accordance with random births and deaths) . Ex- 
amples of substantial alteration of finite systems under ef- 
fect of intrinsic noises gives epidemics |15|- 17 1, predator- 
prey population dynamics [18, 19], opinion dynamics [2o| . 
biochemical clocks J2l|, , genetic networks 23 1 , cyclic 
trapping reactions j24|, etc. Within the phase-plane lan- 
guage, phase transitions pointed out present the simplest 
case, where a fixed point appears only. We are interested 
in studying more complicated situation, when the system 
under consideration may display oscillatory behavior re- 
lated to the limit cycle appearing as result of the Hopf 
bifurcation 25, [2^. It has long been conjectured [23] 
that in some situations the influence of noise would be 
sufficient to produce cyclic behavior [2^. Recent con- 
sideration [29| allows the relation between the stochastic 
oscillations in the fixed point phase and the oscillations 
in the limit cycle phase to be elucidated. In the last case, 
making use of co-moving frame allows fiuctuations trans- 
verse and longitudinal with respect to the limit cycle to 
be effectively decoupled. It appears while the latter fluc- 
tuations are of a diffusive nature, the former ones follow a 
stochastic path. To formulate related model we consider 
the system with a finite number of constituents related 



to components of the state vector [30] 

n = iVX + \/7Vx. (1) 

Characteristically, a deterministic component X is pro- 
portional to total system size N < 00, whereas a random 
one X is the same to its square root. In the limit of in- 
finite particle numbers N ^ 00, such systems are faith- 
fully described by deterministic equations to find time 
dependence X(t), which addresses the behavior of the 
system on a mean-field level. On the other hand, a sys- 
tematic study of corrections due to finite system size can 
capture the behavior of fluctuations x(t) about the mean- 
field solution. These fluctuations are governed with the 
Langevin equations, however, in difference of approach 
[2^ , we consider multiplicative noises instead of additive 
ones, on the one hand, and nonlinear forces instead of 
linear ones, on the other. Within such framework, the 
aim of the present paper is to extend analytical descrip- 
tions [2^ of finite-size stochastic effects to systems where 
noises play a crucial role with respect to periodic limit 
cycle solution creation or its supression. We will show 
that character of the system behavior is determined with 
relation between scales of temporal variation of principle 
variables and their coupling. 

The paper is organized along the following lines. In 
Section |lll we obtain conditions of the limit cycle cre- 
ation making use the pair of stochastic equations with 
nonlinear forces and multiplicative noises. Sections IIIH 
llVl are devoted to consideration of these conditions on the 
basis of stochastic Lorenz system with different regimes 
of principle variables slaving. According to Section IIIII 
the limit cycle is created only in the case if the most fast 
variation displays a principle variable, which is coupled 
nonlinearly with two degrees of freedom or more. Oppo- 
site case is studied in Section IIVI to show that the limit 
cycle disappears in non equilibrium steady state. Section 
IVl concludes our consideration. 
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II. STATISTICAL PICTURE OF LIMIT CYCLE 

According to the theorem of central manifold ^25], to 
achieve a closed description of a limit cycle it is enough 
to use only two variables Xa, a — 1,2. In such a case, 
stochastic evolution of the system under investigation is 
defined by the Langevin equations [s^ 

ia = /("^+SaCW, a =1,2 (2) 

with forces /*^"-' — (a^i, a;2) and noise amplitudes 
Ga — Ga{xi, X2), being functions of stochastic variables 
— 1,2, and white noise ((t) determined as usu- 
ally: (C(t)) = 0, (C(OC(^')> = - t')- Our principle 
assumptions are as follows: (i) white noise is equal 
for both degrees of freedom Xa ; (ii) microscopic transfer 
rates are non correlated for different variables Xa ■ Then, 
the probability distribution function V = V{xi, X2]t) is 
determined by the Fokker-Planck equation 

or dJ° 
dt dxa 

where sum over repeated Greek indexes a = 1, 2 is meant 
and components of the probability current take the form 

2 oxfj 



= 0, 



(3) 



with the generalized forces 



dxp 



(5) 



being determined with choice of the calculus parameter 
A S [0,1]. Within the steady state, the components of 
the probability current take constant values Jq"'' and the 
system behaviour is defined by the following equations: 



0x2 



(1) 



(6) 

2jr(2)-p ^ _2J($^'. (7) 



Multiplying the first of these equations by factor Q2 
and the second one by Qi and then subtracting results, we 
arrive at the explicit form of the probability distribution 
function as follows: 



V (xi,X2) = 



J^^Q2 {XI,X2) " Jq^'Qi {xi,X2) 
V{xi,X2) 

V{XI,X2)= (^2^^^^ 

t;^ ) - Q1Q2 ( 
0x2 I \ oxi 



(2). 



dxi 



^1.^(2)) (8) 
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This function diverges at condition 



2 (^1.^(2) - ^2.^(1- 



9^2 ^2 8G1 
^2^ — 

, 0x1 0x2 



^g.g2r—-'-^ 

OX I 0x2 



(9) 



that physically means appearance of a domain of forbid- 
den values of stochastic variables x^, which is bonded 
with a closed line of the limit cycle. Characteristically, 
such a line appears if the denominator 'D{xi,X2) of frac- 
tion ([5]) includes even powers of both variables xi and 

X2. 

It is worth to note the analytical expression ([5]) of the 
probability distribution function becomes possible due 
to the special form of the probability current ^ , where 
effective diffusion coefficient takes the multiplicative form 
Pag = GaGfS- In general case, this coefficient is known 
[33| to be defined with the expression 



(10) 



where kernel lat determines transfer rate between mi- 
croscopic states a and b, whereas factors 5° and are 
specific noise amplitudes of values Xa related to these 
states. We have considered above the simplest case, when 
the transfer rate lab = I is constant for all microscopic 
states. As a result, the diffusion coefficient pUj) takes the 
needed form = QaQp with cumulative noise ampli- 
tudes = VlT,a9a and gf3 = y/lT,b9p- 



III. NOISE INDUCED RESONANCE WITHIN 
LORENZ SYSTEM 

As the simplest and most popular example of the self- 
organization induced by the Hopf bifurcation, we con- 
sider modulation regime of spontaneous laser radiation, 
whose behaviour is presented in terms of the radiation 
strength E, the matter polarization P and the difference 
of level populations S [30| . With accounting for stochas- 
tic sources related, the self-organization process of this 
system is described by the Lorenz equations 



teE = [-E + geP - '^{E)] + gEdt), 
TpP = {-P + apES) + gpCit), 
TsS=[iSe~S)~asEP]+gsCit). 



(11) 



Here, overdot denotes differentiation over time t] te,p,s 
and aE,p,s > are time scales and feedback constants of 
related variables, respectively; gE,p,s are corresponding 
noise amplitudes, and is driven force. In the absence 
of noises {gE = gp = gs = 0) and at relation tp,ts ^ te 
between time scales, the system addresses to limit cycle 
only in the presence of the nonlinear force [s^ l 



kE 



l + E^/El 



(12) 



characterized with parameters k > and i?„. In this 
Section, we consider noise effect in the case of opposite 
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relation te <C tp,ts of time scales, when periodic varia- 
tion of stochastic variables becomes possible even at sup- 
pression of the force 

It is convenient further to pass to dimensionless vari- 
ables t, C, E, P, S, qe, gp, gs with making use of the 
related scales: 



-1/2. 



Ps = (alapas)-^/^. 
5b = (Tp/apas)^^^, gp = 



Es = iapas)-^/^ 
Ss = {aEapy^; 
-- {Tpla\apasf^ , 



(13) 



1/2, 



9s = V ' lasap. 

Then, equations (llip (these equations are reduced to ini- 
tial form 31 1 if one set there X = y^a/eE, Y = ^Ja/eP, 
Z = Se — S, r = Se, and b = cr/e.) take the simple form 

a-^E = -E + P- ^(E) + gECit), 

P=-P + ES + gpat), (14) 

{e/a)S={Se-S)-EP + gsat), 

where the time scale ratios 



tp/te, 



ts/te 



(15) 



are introduced. In the absence of noises, the Lorenz sys- 
tem (jl4l) is known to show the usual bifurcation in the 
point Se — 1 and the Hopf bifurcation at the driven force 



TP r^^ 



TE 



3t, 



(16) 

'''s '^P 

However, at gE ~ gp ^ gs ~ the limit cycle is unstable 
and the Hopf bifurcation arrives at the strange attractor 
only. 

With switching on noises, the condition te ^ Tp al- 
lows for to put l.h.s. of the first equation (fTi|) to be equal 
zero. Then, the strength is expressed with the equality 



E = P + gECit), 



(17) 



whose insertion into the system (ITil) reduces it into two- 
dimensional form 

P = -P{l-S) + GpC{t), 
S = (a/e) [{Se -S)- + GsCit) (18) 
with the effective amplitudes of multiplicative noises 



Gp^^gl + glS^, gs^{Tp/Ts)^gl + glP^ (19) 
and the generalized forces 



-P(l-5) + A 



9l 



■S, 



'{gs/gEy + P^ 



ts/tp y {gp/gEY + S^ ' 
.F(^) - (Tp/rs) [{Se -S)- P2] (20) 



+A 



\gp/gEr + s^ 
ts/tp^ V i9s/gEy + P^' 



9l 



■P 




FIG. 1: Steady-state distribution function ([8]) at Jq 



J 



(S) 



10, TP = Ts, S'e = 0.5, qe = 0.5, gp = 1.376, gs = 2.5 



In this way, the probability density 
ues at condition 



takes infinite val- 




^ + S^P{1-S) 
9e 



9p. + sAM+P^[{Se-S)-P^] (21) 

9 E / \ 9 E 



9%^ f9l p2 
+ 2 



r. 9% ( g'p q2 



p = o. 



where we choose the simplest case of Ito calculus (A = 0). 

Reduced Lorenz system has two-dimensional form 
(|2|), where the role of variables xi and X2 play the mat- 
ter polarization P and the difference of level populations 
S. According to the distribution function ([8]) shown in 
Figdl the stochastic variables P and 5 are realized with 
non-zero probabilities out off the limit cycle only, whereas 
in its interior the domain of forbidden values P, S ap- 
pears. That is principle difference from the deterministic 
limit cycle, which bounds a domain of unstable values of 
related variables. The form of this domain is shown in 
Figl2]at different values of the noise amplitudes g^;, gp, 
gs and driven force Se- It is seen, this domain grows 
with increase of the driven force Se, whereas increase of 
the force ffuctuations gE shrinks it. On the other hand, 
phase diagrams depicted in Figl3]show the strengthening 
noises of both polarization and difference of level popula- 
tions enlarges domain of the limit cycle creation (in this 
way, force noise gE shrinks this domain from both above 
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40 S 



FIG. 2: Form of limit cycle at: a) qe = 0.5, gp — 11, = 6 
(curves 1-3 relate to Se = 0.5,1.0,2.0, correspondingly); 
b) Se = 0.5, gp — 7.5, gs ~ 6.5 (curves 1-3 relate to gE ~ 
1.0,0.6,0.5, correspondingly) 
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FIG. 4: Steady-state probability distribution as function of 
the radiation strength E and the difference of level popula- 
tions S at conditions tp <^ te = ts for parameters k = 10, 
Se = 11.6, gE = 0.2, gp = 0.2, gs = 0.2 and different proba- 
bility currents Jo^\ Jq'^' (shown on the related panels) 



FIG. 3: Phase diagrams of the limit cycle creation deter- 
mined with Ea. (|21|l at: a) gE = 0.5, curves 1-4 correspond 
to Se = 0.0,0.5, 1.0,2.0, respectively; b) Se = 0.5, curves 1-3 
correspond to gE = 0.1, 0.5, 1.0, respectively (diamonds relate 
to the values gp and gs, for which limit cycles in Fig[2]are 
depicted) 



and below, whereas increase of driven force Sf, makes the 
same from above only). 



IV. LORENZ SYSTEM WITHOUT LIMIT 
CYCLE 

At condition rp ^ te,ts, the deterministic system 
{gE,p,s — 0) has a hmit cycle only at large intensity k of 
non linear force p2)) 34| . In this case, it is convenient to 
measure the time t in the scale te and replace rp by te 
in set of scales (IT3l) . Then, one obtains the relation (cf. 
Eq.Cll)) 



P^ES + gpC{t) 



(22) 



and the Lorenz system ([T^ is reduced to two-dimensional 
form 



E = - [E{1 ~S) + ^{E)] + GECit), 
5 = e-i [Se~S{l + E^)] +gsC{t) 

with the effective noise amplitudes 

Ge-^ 



gl+9% gs = e-^Jgl+glE\ 



(23) 



(24) 



The generalized forces are as follows: 

.F(^) = e-^ [{S, -S)- SE^] 



(25) 



^X^E 



1 + {gE/gp? 



e \l [gslgpY + E^- 

The probability distribution function ([5]) diverges at con- 
dition 

(^^/^^^l^[^(i?)+i?(l-5)] 



1 + {gE/gpY 
l igs/gpf + E'^ 
1 + {gE/gpf 

+ ( A 



[S, - S{1 + E^)] 

^ ' g'pE^ 0, 



(26) 



being the equation, which does not include even powers 
of the variable S. 

As a result, one can conclude that departure from equi- 
librium steady state destroys deterministic limit cycle at 
the relation rp ^ te,ts between characteristic scales. 
This conclusion is confirmed with Fig|3] that shows di- 
vergence of the probability distribution function on the 
limit cycle of variation of the radiation strength E and 
the difference of level populations S at zeros probability 
currents Jq^"* and Jq^^ only. With increase of these cur- 
rents the system escapes from equilibrium steady state 
and maximum of the distribution function shifts to non- 
closed curves being determined with equation 
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V. CONCLUSION 

We have considered effect of stochastic sources upon 
self-organization process being initiated with creation of 
the Hmit cycle induced by the Hopf bifurcation. In Sec- 
tions IllllllVi we have applied general relations obtained 
in Section to the stochastic Lorenz system to show 
that departure from equilibrium steady state can destroy 
or create the limit cycle in dependence of relation be- 
tween characteristic scales of temporal variation of prin- 
ciple variables. 

Investigation of the Lorenz system with different re- 
gimes of principle variables slaving shows that additive 
noises can take multiplicative character if one of these 
noises has many fewer time scale than others. In such a 
case, the limit cycle may be created if the most fast vari- 
able is coupled with more than two slow ones. However, 



the case considered in Section |IV] shows such dependence 
is not necessarily to arrive at limit cycle, as within adi- 
abatic condition Tp <^ te, both noise amplitude Gs{E) 
and generalized force T'^^\E), determined with Eqs. 
(IM)) . (pS)) . are functions of the squared strength only. 
The limit cycle is created if the fastest variations dis- 
plays a principle variable, which is coupled with two dif- 
ferent degrees of freedom or more. Indeed, at the relation 
te Tp, Ts of relaxation times considered in Section [TTTl 
variations of the strength E in nonlinear terms of two 
last equations arrive at double- valued dependencies 
of the noise amplitudes Gp and Gs on both difference of 
level populations and polarization, which are appeared 
in Eqs. (fT9|) as squares and P^. This appears physi- 
cally as noise induced resonance related to the limit cycle 
created by the Hopf bifurcation, that has been observed 
both numerically 13[ and analytically 35 1. 
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